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Abstract 

We show that, under certain combinations of the parameters governing the interaction of a 
harmonically trapped ion with a laser beam, it is possible to find one or more exact eigenstates of 
the Hamiltonian, with no approximations except the optical rotating-wave approximation. These 
are related via a unitary equivalence to exact eigenstates of the full Jaynes-Cummings model 
(including counter-rotating terms) supplemented by a static driving term. 
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I. INTRODUCTION 



In recent years, trapped ions interacting with laser beams [1, 2] have become an extremely 
interesting system for the investigation of fundamental physics. For instance, they have 
been used to generate nonclassical motional states [3, 4, 5], and have been proposed for 
applications such as quantum computation [1, 6, 7, 8] and precision spectroscopy [9]. 

Despite the relative simplicity of this system, the full theoretical treatment of its dy- 
namics is a nontrivial problem, as the laser-ion interaction is highly nonlinear. Even in 
the simplest case where only a single ion is in the trap, one is usually forced to employ 
physically motivated approximations in order to find a solution. A well-known example is 
the Lamb-Dicke approximation [11], in which the ion is considered to be confined within 
a region much smaller than the laser wavelength. Many treatments also assume a weak 
coupling approximation, i.e., a sufficiently weak laser-ion coupling constant. Under these 
conditions, tuning the laser frequency to integer multiples of the trap frequency results in 
effective Hamiltonians of the Jaynes-Cummings type [12], in which the centre-of-mass of the 
trapped ion plays the role of the field mode in cavity QED [13, 14]. 

Recently, a new approach to this problem has been suggested [10], based on the applica- 
tion of a unitary transformation T which linearises the total ion-laser Hamiltonian. Under 
this transformation the Hamiltonian becomes exactly (not effectively!) equivalent to the full 
Jaynes-Cummings model (JCM), inchiding counter- rotating terms, together with an extra 
atomic driving term. Remarkably, the ion-trap system is thus formally equivalent to an 
atom interacting with a single-mode quantised electromagnetic field. 

The existence of this correspondence is very useful, since it allows us to map interesting 
properties of each model onto their counterparts in the other. For example, we can make use 
of the well-known 'rotating- wave' limit where the JCM is analytically, albeit approximately, 
soluble [12]. Using the map T"^ to translate this solution back into the ion-trap scenario 
has led to the identification of a new dynamically interesting regime for that system [10], 
where phenomena such as "super-revivals" in the ion-laser interaction can occur. The same 
solution has also led to a scheme for realising substantially faster logic gates for quantum 
information processing in a linear ion chain [8]. 

In the present note, we apply this line of reasoning to the problem of obtaining exact 
eigenstates of either system. To our best knowledge, this has never been achieved in the 
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case of the ion-trap Hamiltonian. In the case of the JCM (with no extra driving), exact 
eigenstates have been derived for those combinations of system parameters at which level 
crossings occur in the spectrum [15, 16]. In this paper we use a simple physical ansatz similar 
to the one in [16] to derive a family of exact eigenstates for the ion-trap system (again, these 
are valid only under certain combinations of system parameters). In the special case where 
the ion is driven on resonance, these states correspond, via the mapping T, to the ones found 
in refs. [15, 16]. In addition, we study the asymptotic limit of large Lamb-Dicke parameter rj 
in the ion-trap system, which under T corresponds to a large coupling constant in the JCM. 
We obtain in this way a very simple explanation of the asymptotic forms of all eigenvalues 
and eigenvectors in this limit (these have also been derived in [16] using other methods). 

II. ION TRAP-JCM CORRESPONDENCE 

Let us start then by recalling the equivalence between the ion-trap system and the JCM 
[10]. The Hamiltonian for the ion-laser dipole interaction, with no approximations (except 
the optical rotating- wave approximation [9]) can be written as 



5 = uJatom~^iaser the laser-iou detuning, fl the (real) Rabi frequency of the ion-laser coupling, 
1] the Lamb-Dicke parameter, and we have chosen units where h = 1. On the other hand, 
the Jaynes-Cummings Hamiltonian with counter-rotating terms can be written as 



Although these two models appear to be physically and mathematically quite distinct, 
they are in fact exactly equivalent The easiest way to see this is by rewriting eq.(l) in a 
notation where operators acting on the internal ionic levels are represented explicitly in 
terms of their matrix elements: 





(2) 




(3) 
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Consider now the unitary operator 

^i(n'W 6(0) 

V2y-Dl(l3) D(I3) 
where /5 = ir]/2. It is possible to check after some algebra [10] that: 

Returning to the usual notation, we obtain 

7i = i/n + 0(7;, + ^ (cr+ + (7_) (a - a^) + ^ {a+ + (7_) + ^ (6) 

Comparing with eq. (2) it can be seen that this is precisely the Jaynes-Cummings in- 
teraction, supplemented by two additional terms: the first corresponds to an extra static 
electric field interacting with the atomic dipole, and the second is just a constant energy 
shift. In particular, a purely Jaynes-Cummings form is recovered when S — 0, correspond- 
ing to a resonant laser- ion interaction in eq. (1). Of course, the various parameters of the 
Hamiltonian in eq. (6) have different meanings than they do in eq. (1): u becomes the cavity 
field frequency uj, 2Q the atomic transition frequency ujq, 6 the coupling strength with the 
static field and 77 the ratio between the Jaynes-Cummings Rabi frequency 2A and the cavity 
frequency cu. In what follows we shall refer to eq. (6) as the 'Jaynes-Cummings picture' of 
the ion-trap Hamiltonian, eq. (1). 

As noted above, this correspondence enables one to map interesting properties of each 
model onto the other. In this paper, we will use it to translate between corresponding 
eigenstates of each system (it is clear that, if {ip) is an eigenstate of Hion, then T {ip) is a 
corresponding one for H.). In this regard it is important to point out that, although Hion and 
H both describe systems consisting of a two-level atom interacting with a bosonic mode, one 
should not identify each of these subsystems with their counterparts after the transformation 
has been applied. This is due to the fact that T is an entangling transformation: separable 
internal-motional states of the trapped ion can be mapped into entangled atom-cavity states 
in the corresponding cavity QED system. 
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III. EIGENSTATES IN THE ASYMPTOTIC LIMIT 



In order to demonstrate the convenience of this correspondence, let us consider what 
happens to the systems above in the 'asymptotic' hmit where 77 ^ Jl/i/ (this problem has 
been studied in [16] in the special case 5 — 0, using different methods). In the ion-trap 
picture, this limit corresponds to having a very large atomic recoil after the absorption of 
a photon from the laser; in the JCM picture, to having a very large atom-mode coupling 
constant. From the point of view of current ion-trap and cavity QED experiments, neither of 
these conditions can be considered realistic. Nevertheless, it may well be possible to engineer 
such a situation artificially, for example an 'effective' JCM with large 'couphng constant' 
can be obtained when a single ion is simultaneously illuminated with weak lasers tuned to 
the first red and blue sidebands [1] . 

It is simpler to consider the problem first in the JCM picture. If 1] ^ Q/u, then the 
term containing fl in eq. (6) makes a negligible contribution to the total energy. We can 
therefore expect that the asymptotic eigenvalues and eigenvectors will not be affected if 
this term is in fact absent. Setting Q = implies that = (t+ + (7_ becomes a constant 
of the motion, so the asymptotic eigenstates can be chosen to be of the form |±) (8) |0±), 
where |±) = (l^f) ± |e)) /\/2. In order to determine their form explicitly, and also their 
corresponding eigenvalues, let us switch back to the ion-trap picture. Setting f2 = in eq. 
(1) reduces the Hamiltonian to that of an uncoupled spin-boson system. The asymptotic 
eigenstates of the ion-trap are therefore simply |em),|g'm). Transforming back to the 
JCM picture using eq. (4), we obtain thus = \±ir]/2;n), where \a;k) = D (a) \k) is 
a displaced number state [18] (note that they are independent of 5). The corresponding 
eigenvalues are mu ±6/2 [21]. Note in particular that, when \S\ = kv for integer k, all but 
k of these energy levels arc degenerate. In the special case 5 = it can in fact be shown 
from symmetry considerations alone that the entire asymptotic spectrum must be two-fold 
degenerate [19]). Conversely, no such asymptotic degeneracies occur for \S\ ^ ku. 

IV. EIGENSTATES FROM A SIMPLE ANSATZ 

Let us return now to the ion-trap Hamiltonian, eq. (1). We will construct an ansate which 
allows the determination of exact eigenstates of this system, provided certain relations are 
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satisfied between the parameters Q, 5, 77. Consider the possibihty of finding an eigenstate of 
the form 



n 



m+l 



n=0 



Let us see whether the eigenvalue equation 

can be satisfied. Eq. (3) shows that this requires |0) to be of the form 



m+l 



m+l 



{ir]) ^ dn \n) = ^ (i„ \-i7]; 



n) 



n=0 



n=0 



We thus require 

m+l 



Hicm 



(7) 



(8) 



(9) 



(("""^ ^) ^""^^"^ I"") + (^-^Cn + (in (^-T,] ]\9) \-W,n). (10) 



Now, using the simple fact that D'^ {a) dD (a) = d + a [20], it is easy to show that 
displaced number states satisfy the recursion relation 



n \a: 



]k) = {\a\^ + k) \a]k) + a\/k + 1 \a]k -\- 1) + a*\fk |q;; A; — 1) . 



(11) 



Substituting then eqs. (7), (10) and (11) into cq. (8) and using the fact that {|q;; A;)}^o 
an orthonormal basis gives the following eigenstate conditions: 

drt. 



m + 1 — n 



; < n < m 



where the dn coefficients satisfy lim+i = and 



IT] Sm-i -irjV2 
ir]\/2 



y/'m + 1 dm', n — 171 + 1 



do 



El —iTj^/m 



d„ 



(12) 



(13) 
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where 

., = (, + l-,^) + _-^_^. (14) 

Using operator T we can map this state into a corresponding eigenstate of the generahsed 
JCM model in eq. (6) 

\^7^m) = T |C:> ^J2[dn\+)- -cn I-) |-W2;n) . (15) 

Once again we verify that the solution for this system can be easily expressed in terms of a 
displaced number state basis. 

Further solutions can also be found if we note that the Hamiltonian Hion is invariant 
under the combined transformations 

|e) ^ \g); 5^-6; rj ^ -rj. (16) 

Defining then 

„ m+l m 

ICn > = - E < 1^) 1^) + \e) H n) , (17) 



n=Q n=0 



and applying this symmetry transformation to eqs. (10)- (14) we can see that [tp"^ ) is also 
an eigenstate of Hicm, with eigenvalue E~ = (m + l)z/ — |, as long as 



< n < m 



m + 1 — n 



(18) 



2 

-\/m-\-l d'^; n — m + 1 



02 

where the d'^ coefficients satisfy an equation analogous to eq.(13) but with 

Condition (13) means that each of these ansatze succeeds only for certain combinations 
of fl, 6, Tj. This is because the tridiagonal matrix above (which we will refer to as M^) must 
have a zero eigenvalue, or equivalently det — 0. Since det is a polynomial of degree 
m + 1 in ^2 (5 or 772, fixing two of these quantities determines up to m + 1 real solutions for 
the third one. For m < 3 it is possible to solve this 'compatibility' condition algebraically, 
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resulting in explicit expressions in terms of any of the three parameters. For example, in 
terms of 77, the first few solutions are 

m = 0: T]^ ^1-a + sd (20) 
m = 77^ = 2 + sd- — ±W— -- + 2 + sd (21) 

where a — {ft/v)'^ , d — S /u and s = ±1 according to which ansatz we are referring to. For 
m > 3 a numerical solution is necessary. 

The physical meaning of these solutions can be clarified by examining the spectrum of 
Hion, obtained by numerically diagonalising this operator (or, equivalently, Ti). Let us first 
examine the case where 5 = (previously studied in [15, 16]). In fig. (la) we plot the first 
few energy levels as functions of r], for the case where fl — i//2 (corresponding to a resonant 
JCM with no extra driving in eq.(6)). It can be readily seen that pairs of adjacent energy 
levels form braids around the lines E — mu for integer m, converging to these lines in the 
limit of large 1] (as expected from our considerations in section III). The level crossings 
occur precisely on top of these lines, and there are m such degeneracy points on the line 
E = mz/ [21]. (In fact, it can be shown [16] that these are the only energy values at which 
degeneracies may occur in this system). 

Note now that, when 5 — equations (14) and (19) coincide, and thus, 
simultaneous (in fact, degenerate) eigenstates for any combination of Q, 77 at which det — 
0. Another way of understanding this fact is to note that the ion-trap Hamiltonian Hion has 
in this case an extra symmetry: it commutes with the parity-like observable exp^ina'' a) 
(in the JCM picture, the corresponding symmetry operator is cxp[i7ra^ a)) . It can be easily 
checked that neither IV'Sn ) nor \^]2n) have this symmetry, but simple linear combinations 
of them do. Every solution to the equation det — must therefore correspond to one 
of the level crossing points; in other words, the location of these points can be calculated 
directly from this condition (for m > 3 this cannot be done algebraically, of course). In fact, 
it turns out that, when Q = z//2, this equation always has m real and positive solutions for 
1] [16], and so every crossing point in fig. la is accounted for in this way. (See [16] for a 
detailed discussion of what happens as the value of fl is changed) . 

For more general values of S the situation changes as follows (fig lb,c): whenever S/u — k 
where A; is a nonzero integer (i.e., the laser is tuned to a sideband transition), it is easy to see 
that E^ — E^^f^ = (m+l+k/2)iy. In other words IV'I^J') and ~) are again degenerate 
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eigenstates, and again the degenerate energies coincide with the asymptotic values for large 
77 (recall sec. III). In fact, it can be verified that detM+ = =^ detM~_^j. = 0, so that 
both eigenstates occur for the same system parameters, thus corresponding once again to 
line crossings in the spectrum (fig. lb). (For low values of m, k this can be shown explicitly, 
for higher values one can again resort to numerical methods). It seems likely that this 
coincidence may once again be due to an underlying symmetry, however in this case we 
have not been able to determine it. Finally, when 5 is not an integer, the crossings in the 
spectrum become avoided crossings (fig. Ic), and in fact it is easy to show, using a method 
analogous to that in ref. [16], that no degeneracies can occur in the system. In this case, 
the solutions obtained above simply mark (some of) the points where the spectrum crosses 
its asymptotic value (fig. Ic). 



V. CONCLUSION 



In conclusion, we have shown that in certain circumstances it is possible to obtain exact 
eigenstates for the Hamiltonian of a single trapped ion. These states are often (but not 
always) connected with the existence of level crossings in the system's energy spectrum 
viewed as a function of the Lamb-Dicke parameter. Another property is that they are always 
expressible as a finite expansion in terms of certain physically well-motivated vibrational 
states (number states and displaced number states). In this regard, it is worth recalling 
a long-standing conjecture concerning a possible complete solution to the JCM in terms of 
known functions (valid for any values of the system parameters). It was suggested in [22] 
that the eigenstates of this system may always be expressible as a finite expansion in terms 
of a certain basis set, which in Bargmann representation corresponds to a particular class 
of transcendental functions. As far as we know this proposal has never been analytically 
verified, although numerically it seems to hold [22, 23]. In the light of our present results, 
it is natural to speculate that this result should still hold in the presence of an additional 
static field, and thus should also apply to the single-ion system. 

H. M.-C. would like to thank CONACYT (Consejo Nacional de Ciencia y Tecnologia) and 
the Royal Society of London for support, and also Imperial College for hospitality. D.J. was 
supported in part by the Brazilian Government's Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico (CNPq) and by Pembroke College, Cambridge. 
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FIG. 1: First few energy levels (in units of h) of a single trapped ion driven by a laser beam, as 
a function of the Lamb-Dicke parameter rj, for coupling constant 0, = v jl. Dark circles denote 
parameter values for which exact eigenstates of this system can be found using the ansatz described 
in section IV. (a) When the laser is resonant ((5 = 0), a solution exists whenever the energy E is 
an integer multiple of the trap frequency z^, or in other words whenever the curves intersect the 
horizontal asymptotes to which all energy levels tend in the limit 77 — > oo. These arc also the level 
crossing points in the spectrum [16]. (b) A similar behaviour occurs when 5/u = k for integer 
A; 7^ (i.e., the laser is tuned to one of its sidebands) (here k = 1). In this case the level crossings 
are at points where E = {m + 1 + 0.5)i/, which are again the asymptotic values of E. (c) Finally, 
when S is not tuned to a sideband {S = O.bu here), level crossings become avoided crossings, and 
our ansatz finds exact eigenstates only at some of the points where the curves cross the asymptotes 
Em = {m + 1.25)u. 
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